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^ ! ABSTRACT 

o 

Active galactic nuclei (AGNs) are characterized by a clear correlation between luminosity and metallicity (£agn _ Zagn relation). 
The origin of this correlation is not clear. It may result from a relation between the black hole mass (M B h) and metallicity, or 
\^») , from a relation between the accretion rate (L/L Edd ) and metallicity. To investigate the origin of the Lagn-Zagn relation, we use 

■ optical spectra of 2383 quasars at 2.3 < z < 3.0 from the Sloan Digital Sky Survey. By using this data set, we have constructed 
composite spectra of 33 subsamples in intervals of both M BH and L/L E(]d . From these composite spectra we measure emission-line 
flux ratios that are sensitive to the metallicity of the broad line region (BLR); specifically, N V/U240/C iv/11549, N V/U240/Heii/11640, 
(SiiV/tl398+Oiv]/U402)/CiV/tl549, and Alni/I1857/C ivil549. We find that there is a significant correlation between M BH andZ BLR 
as inferred from all four metallicity-sensitive emission-line flux ratios. This result strongly suggests that the observed Lagn-Zagn 
relation is mostly a consequence of the M bh -Z A gn relation. The relation between M BH and Z BLR is likely a consequence of both the 
M B H-Mb„i relation and of the mass-metallicity relation in the host galaxy. We also find that L/Leu correlates with the emission line 
flux ratios involving N v (more specifically, N v/C iv and N v/He n), while it does not correlate with the other two metallicity sensitive 
emission line flux ratios, i.e., (Siiv+O iv])/C iv and Al m/C iv. These correlations indicate that the emission-line flux ratios involving 
N v depend on both metallicity and relative abundance of nitrogen. We suggest that the relation between L/L Ei(i and those line ratios 
involving nitrogen, is caused by a delay of the black hole accretion rate relative to the onset of nuclear star formation of about 10 s 

C/3 , years, which is the timescale required for the nitrogen enrichment. 

a , 
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■ 1. Introduction past star forming history of the host galaxy (e.g., Matteucci & 
\ Padovani 1993; Hamann & Ferland 1993, 1999; Venkatesan et 

(X) , In the local universe, there is a tight correlation between the al. 2004). Therefore, it is possible to investigate the formation 

• mass of black holes at the center of galaxies (M B h) and their and evolution of SMBHs and of their host galaxies by compar- 

_J | hosting stellar spheroids (M bu i), as well as with their velocity ing L A gn and the gas metallicity. However it is not clear how 

, dispersion (e.g., Kormendy & Richstone 1995; Richstone et al. the gas metallicity of BLR is connected to the chemical prop- 

(-->>, 1998; Ferrarese & Merritt 2000; Gebhardtet al. 2000; Tremaine erties of their host galaxies. The BLR is a very small region 

y—i ! et al. 2002; Marconi & Hunt 2003; Haring & Rix 2004; Barth et in galactic nuclei (r BL R < 1 pc; e.g., Suganuma et al. 2006), 

^ . al. 2005). These correlations suggest that the formation and evo- which may have evolved more rapidly than their host galaxies, 

lution of supermassive black holes (SMBHs) are connected with However, a clear relation between the metallicity of narrow-line 

those of their host galaxies. However the mechanism responsible region (NLR) and the AGN luminosity has been also reported in 

$J] ■ for such co-evolution is not yet clear. In this study, we approach recent studies (e.g., Nagao et al. 2006b; Matsuoka et al. 2009). 

Ci | this problem by focusing on the properties of the active galactic In contrast to the BLR, the scale size of the NLR is comparable 

nuclei (AGNs), in particular M BH and accretion rate L/L Edd , and to the size of the host galaxies (r NLR ~ 10 2 ~ 4 pc; e.g., Bennert et 

their gas metallicity (Z A gn)- al. 2006a, 2006b) and its gas mass is about M NLR ~ 10 5 ~ 7 M o 

Hamann & Ferland (1993) proposed that emission-line flux ( much lar S er than for the BLR < m blr ~ 1O 2 - 4 M ; see, e.g., 

ratios Nv^l240/Civil549 and NvA1240/He nvll640 trace the Baldwin et al. 2003), implying that the NLR metallicity traces 

broad-line region (BLR) metallicity (Z BL r) across a wide range the enrichment in host galaxy. As a consequence, the finding of 

of the AGN luminosities (L AGN ) and BLR densities. They also a Ugn-Znlr relation suggests that also the L agn -Zblr relation 

reported a clear relation between the AGN luminosity and BLR traces a connection between the nuclear engine and the metallic- 

metallicity: more luminous AGNs tend to have more metal-rich ^ ln the host g alax y- or at least ln the central re 8 lon - 

BLR clouds (see also Nagao et al. 2006a; Juarez et al. 2009). Since in flux-limited samples Mbh and L A gn are correlated, 

The AGN luminosity is determined by both the black hole mass Hamann & Ferland (1993) suggested that the L A gn-Zblr relation 

and accretion rate relative to its Eddington limit. On the other is a consequence of a relation between black hole mass and BLR 

hand, the metallicity is an important tracer for the formation and metallicity. However, the L A gn-Zblr relation may also originate 

evolution of the host galaxy, since it is closely related to the from a connection between the gas metallicity and the accretion 
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rate L/LEdd- Warner et al. (2004) analyzed rest-frame ultravio- 
let spectra of type 1 AGNs at < z < 5 and found that the 
emission-line flux ratios involving nitrogen lines (expected to 
trace the BLR metallicity; see e.g., Hamann et al. 2002) show a 
correlation with the black hole mass, while no correlation with 
the Eddington ratio. On the other hand, Shemmer et al. (2004), 
by using near-infrared spectroscopic measurements of luminous 
quasars at 2 < z < 3.5, found that the BLR metallicity, as in- 
ferred by N v/C iv, is primarily correlated with the Eddington 
ratio and not with the black hole mass. Dietrich et al. (2009) 
studied 10 high-z quasars and showed that the N v/Civ flux ra- 
tios indicate super-solar metallicities in the quasars with a high- 
accretion rate, suggesting a correlation similar to that found by 
Shemmer et al. (2004), although from those quasars alone there 
is no correlation between N v/Civ and L/L E dd (see their Figure 
12). This discrepancy may be attributed to the small size of the 
adopted samples (578 quasars for Warner et al. 2004 and 92 
quasars for Shemmer et al. 2004). 

We use a large quasar sample of the Sloan Digital Sky Survey 
(SDSS; York et al. 2000) to investigate possible correlations be- 
tween metallicity, black hole mass, and Eddington ratio. We ap- 
ply a stacking method delivering composite spectra by dividing 
the quasar sample into several bins of Eddington ratio and of 
black hole mass. This method allows us to measure emission- 
line flux ratios unambiguously by using the high signal-to-noise 
ratio spectra and to break the degeneracy between black hole 
mass and Eddington ratio in flux limited samples. More specifi- 
cally, we present the composite spectra for 33 subsamples binned 
by M B h and L/L E dd out of 2383 quasar spectra obtained from 
the SDSS, and investigate possible dependence of metallicity- 
sensitive emission-line flux ratios on M B h and L/Lem- In sec- 
tion 2, we describe the quasar sample and composite spectra. 
The description of the measurement of emission lines is given in 
section 3. In section 4 we summarize our results, and in section 
5 we discuss the interpretation for the our results. The summary 
of this paper is given in section 6. 

2. Sample and composite spectra 

We select a large sample of SDSS quasars from the catalog in 
Shen et al. (2008) on the basis of the following criteria: 

1. redshift in the range 2.3 < z < 3.0 to investigate rest-frame 
ultraviolet spectra including emission lines such as N v, C rv, 
and He n, 

2. both the luminosity at 1 350A and the line width of C iv emis- 
sion are measured, and 

3. identified as non broad absorption line (non-BAL) quasars. 

Then we obtain a sample that consists of 2678 SDSS quasars. 
Their spectroscopic data are available in the SDSS archive, Data 
Release 5 (Adelman-McCarthy et al. 2007; Schneider et al. 
2007). The spectral resolution of the SDSS spectra is ~ 2000. 
The wavelength coverage is 3800A :£ /l b s ^ 9200A, within 
which we can measure some redshifted ultraviolet emission lines 
(i.e., N vil240, C iv/11549, He nil640, Si iV/I1398+Oiv]/U402, 
and Al m,U857) to study the metallicity of BLRs. 

The black hole mass is estimated by the line width of 
Crv/11549, adopting the following formula (Vestergaard & 
Peterson 2006; Shen et al. 2008): 



log 



'Mr 



Mr 



= 0.66+0.53 log 



AL A (1350A) 



„, /FWHM\ /1X 
+2l ° g km^T Id) 



10 44 ergs s 

where FWHM is the full width at half maximum of the C iv line. 
We also use the continuum luminosity to estimate the Eddington 
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Fig. 1. Distributions of bolometric luminosity, black hole mass, 
and the Eddington ratio for our quasar sample that of 2678 SDSS 
quasars. 

ratio, L/LEdd with the bolometric correction factor at BC = 3.81 
(Richards et al. 2006; Shen et al. 2008), and we obtain the fol- 
lowing formula, 



log = 5.82+0.47 log 

VLfidd, 



AL A (1350A) 



10 44 ergs s" 



„, /FWHM\ ^ 



Figure 1 shows the distribution of the luminosity, black hole 
mass, and Eddington ratio for our SDSS quasar sample. Our 
sample includes the luminosity at 45.9 < logLb i(ergs s _1 ) < 
48.1, black hole mass at 7.7 < log(M BH /M ) < 11.3, and 
Eddington ratio at -2.6 < log(L/LEdd) < 0.8. Their averages and 
standard deviations are 46.90+0.07, 9.34+0. 1 8, and -0.55 +0. 1 1 , 
respectively. Note that the low cutoffs in the Lb i, M B h, and 
L/L E dd distributions presented in Figure 1 are due to the flux- 
limited nature of the SDSS quasar sample. 

To investigate the relationship between the black hole mass 
and metallicity, independently of the Eddington ratio, or the 
L/L E dd-ZBLR relation independently of M BH , we divide our sam- 
ple into subsamples by fixing one of two parameters of the 
black hole mass and Eddington ratio. By using these subsam- 
ples, we examine possible relationship between the black hole 
mass and BLR metallicity with a given Eddington ratio, or the 
L/L E dd-ZBLR relation at a given black hole mass. Table 1 gives 
the detailed information for our subsamples. Note that we ex- 
clude some quasars that are at the ranges of the black hole mass, 
log(M BH /M ) < 8.4 and 10.0 < log(M BH /M ), and Eddington 
ratio, log(L/L E dd) < -1-2 and 0.0 < log(L/L E ddX and also ex- 
clude some subsamples composed with less than 10 quasars in 
order to improve the statistical significance for results. This im- 
plies a final selection of 2383 quasars. 

We extract the composite spectra for the 33 subsamples 
given in Table 1 in the following way. By adopting the red- 
shift defined in the catalog in Shen et al. (2008), we shift the 
quasar spectra to the rest frame before stacking. After shifting, 
each spectrum is normalized by using the mean flux density 
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Table 1. Definition of the 33 intervals of M BH and L/L Edd used to 
create the composite spectra, along with the number of quasars 
within each interval. 



Tk 1 T T T 1 A If 

Black Hole Mass 


Eddington Ratio 


N 


8.4 < log(M BH /M G ) < 8.6 


-0.4 < log(L/L Edd ) 


< -0.2 


19 




-0.2 < log(L/L Edd ) 


< 0.0 


32 


8.6 < log(M BH /M ) < 8.8 


-0.6 < log(L/L Edd ) 


< -0.4 


14 




-0.4 < log(L/L Edd ) 


< -0.2 


44 




-0.2 < log(L/L Edd ) 


< 0.0 


58 


8.8 < log(M BH /M n ) < 9.0 


-0.8 < log(L/L Edd ) 


< -0.6 


16 




-0.6 < log(L/L Edd ) 


< -0.4 


45 




-0.4 < log(L/L Edd ) 


< -0.2 


71 




-0.2 < log(L/L Edd ) 


< 0.0 


54 




-1.0 < log(L/L Edd ) 


< -0.8 






-0.8 < log(L/L Edd ) 


< -0.6 


71 




-0.6 < log(L/L Edd ) 


< -0.4 


87 




-0.4 < log(L/L Edd ) 


< -0.2 


108 




-0.2 < log(L/L Edd ) 


< 0.0 


55 




-1.0 < log(L/L Edd ) 


< -0.8 






-0.8 < log(L/L Edd ) 


< -0.6 


108 




-0.6 < log(L/L Edd ) 


< -0.4 


190 




-0.4 < log(L/L Edd ) 


< -0.2 


118 




-0.2 < log(L/L Edd ) 


< 0.0 


11 




-1.2 < log(L/L Edd ) 


< -1.0 


22 




-1.0 < log(L/L Edd ) 


< -0.8 


81 




-0.8 < log(L/L Edd ) 


< -0.6 


214 




-0.6 < log(L/L Edd ) 


< -0.4 


188 




-0.4 < log(L/L Edd ) 


< -0.2 


46 


9.6 < log(M BH /M Q ) < 9.8 


-1.2<log(L/L Edd ) 


< -1.0 


25 




-1.0<log(L/L Edd ) 


< -0.8 


126 




-0.8 < log(L/L Edd ) 


< -0.6 


196 




-0.6 < log(L/L Edd ) 


< -0.4 


74 




-0.4 < log(L/L Edd ) 


< -0.2 


18 


9.8 < log(M BH /M Q ) < 10.0 


-1.2<log(L/L Edd ) 


< -1.0 


52 




-1.0<log(L/L Edd ) 


< -0.8 


82 




-0.8 < log(L/L Edd ) 


< -0.6 


61 




-0.6 < log(L/L Edd ) 


< -0.4 


29 



Total 



2383 



at 1350A < 

^rest < 1365 A, where there are no strong emis- 
sion lines. For the stacking procedure, we use the IRAF task 
scombine and make median spectra (instead of average spectra) 
for combining pixels with a 3 sigma-clipping rejection. There 
are two techniques to generate composite spectra "in general": 
(1) the median (or mean) spectrum, which preserves the relative 
fluxes of the emission features; and (2) the "geometric" mean 
(i.e., logarithmic mean) spectrum, which preserves the global 
continuum shape (see, Vanden Berk et al. 2001). Since we are 
interested in the emission-line flux ratios, we apply the former 
method (i.e., the linear median spectrum). Figures 2 and 3 show 
the extracted composite spectra for each bin of black hole mass 
or Eddington ratio, respectively. 



3. Emission-line measurement 

It is sometimes difficult to measure broad emission-line fluxes 
because some close emission-lines (e.g., Lya and N vil240) are 
heavily blended and because the continuum level is often not 
easy to estimate due to the contribution of both weak emission 
lines and Fen multiplet features. Previous studies have mostly 
used one of the two following methods. One is to fit the detected 



emission-line features with some appropriate profile function; 
the other method consists in defining a local continuum level for 
each emission line and integrating the flux above the adopted 
continuum level. The former method is far more accurate when 
the flux measurements requires the deblending of emission lines, 
which is the case for most quasar UV spectra. We therefore adopt 
the former method in this study. 

We use the following function to fit the line profiles: 



F) = 



'Mir 



for A > 
for A < Aq, 



(3) 



where the power-law indices (a and /3) are different between 
the blue side and red side of a given emission line (i.e., a + 
(3 generally). This function achieves a better fit than double- 
Gaussian and modified (i.e., asymmetric) Lorentzian methods 
(see Nagao et al. 2006a for more details). Since the emission- 
line profile depends on the ionization degree (e.g., Gaskell 1982; 
Wilkes 1984, 1986; Baskin & Laor 2005), we divide ultravio- 
let emission lines into the following two groups: high-ionization 
lines (HILs; Nvil240, Oivil402, Niv]il486, Civ^l549, 
and He H/11640) and low-ionization lines (LILs; Sinil263, 
Siiv/11398, Om]il663, Al il11671, A1iilU857, SiiLll887, and 
Cm]/ll909), with the boundary ionization potential of ~ 40 eV 
(e.g., Collin-Souffrin & Lasota 1988 for discussion on this di- 
chotomy). We assume that all emission lines in the same group 
have the same power-law indices of a and /?, and the same wave- 
length shift of the emission line center {Aq) from the systemic ve- 
locity. The line profiles for the same group are fit simultaneously 
to infer the best values of a, f3, and Aq. The absorption by inter- 
vening intergalactic matter on the blue side of Lya affects the 
profile fitting for Lya. Therefore, for the Lya fitting, we adopt 
the red-side index a inferred for HILs, and the other two param- 
eters (J3 and Aq) are left to be free. The model fitting is performed 
in the spectral region at 1210A < i rest < 2000A. The slope and 
the amplitude of the power-law continuum are initially estimated 
in two wavelength regions where the emission line contribution 
to the total flux appears to be small (1445 A < A KSt < 1455 A and 
1973A < Arest < 1983A) and finally determined by the model 
fitting. The fitting results are shown in Figure 4. The best fit 
power-law indices, a and and the systemic velocity of each 
line are given in Table 2. 

The Oi+Sin composite at 1305 A and Cii/11335 are mea- 
sured simply by summing up all of fluxes above the continuum 
level for each line (1286 A < A rest < 1322A for the O i+Si n com- 
posite and 1322A < A res t < 1357 A), because their velocity pro- 
files are different from both HILs and LILs. In each spectrum we 
detect a "1600A bump". This unidentified emission feature has 
been noted in earlier studies, e.g., Wilkes (1984), Boyle (1990), 
Laor et al. (1994), and Nagao et al. (2006a). The possible nature 
of this feature has been discussed in Nagao et al. (2006a) and it 
will not be discussed further in this paper. The spectral regions 
discussed above (Oi+Sinil305, Cii/11335, and 1600A bump) 
are excluded from our fitting procedure. The wavelength region 
at 1687A 

< ^rest < 1833 A is also excluded because it includes 
heavily blended emission lines such as NiV/11719, Alii/11722, 
Nm]/11750, and Fen multiplets, which are very difficult to de- 
blend. Nin]/I1750 would also be useful to estimate the metal- 
licity in addition to Nvil240 (Shields et al. 1976; Baldwin & 
Netzer 1978; Hamann et al. 2002; Dietrich et al. 2003), but we 
did not take it into account because of the difficulty in deblend- 
ing it from broad Fe n emission. 
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-1.2Slog(L/L Edd )<-1.0 




aaog(U m /M e )<10.0 
6Slog(Jf BH /Jf )<9.8 



4Slog(tf BH /tf )<9.6 



-0.8£log(L/Z, Edd )<-0.6 



0.4glog(L/L Edd )<-Q.g 




9.6Slog(Jf BH /Jf )<9.8 
9.4Slog(tf BH /tf )<9.6 



9.2£log(tf BH /J/ G )<9.4 



9.0Slog(y BH /Af o )<9.S 
B.8aog(lf m /Jf a )<9.0 



B.6Slog(Jf BH /Jf )<8.8 



B.4Slog(Jf BH /Jf )<8.6 



rest 



(A) 



-1.0SIog(i/i Edd )<-0.B 




g.ssiog^/jfj.xio.o 

9.6Slog(Jf BH /Jf )<9.8 
9.4Slog(tf BH /tf )<9.6 

__yV_9^i°8(^ B „/Jf G )< 9 - 4 
g.osiogf^/^xg.a 



0.6glog(Z,/L Edd )<-0.4 




9.8Slog(i/ B „/Jf e )<10.0 
9.6Slog(Jf BH /J/ )<9.8 
9.4Slog(tf BH /tf )<9.6 
9.2£log(tf BH /l/ G )<9.4 



9.0 alog(y BH /if o )<9.a 
B.SSlogfAfjn/lfJO.O 



B.6Slog(Jf BH /Jf )<8.8 



0.aglog(L/L Edd )<0.0 




9.2£log(tf BH /l/ G )<9.4 



yNajSloBtJfaa/JfsXQ.O 
/V_J^°s(W"o)<8-B 



B.4Slog(Jf BH /Jf )<8.6 



rest 



(A) 



Fig. 2. Each panel shows the composite spectra in a fixed interval of Eddington ratio by varying Mbh, normalized by the mean flux 
density at 1350A < A KSt < 1365 A. The height scale is the same for all spectra, but different spectra are shifted by two. 



4. Results 

Based on our fitting results described in section 3, we derive 
the emission-line flux ratios that are sensitive to the BLR metal- 
licity, i.e., Nv/Civ, Nv/Hen, (Siiv+Oiv])/Civ, and Alra/Civ 
(Hamann et al. 2002; Dietrich et al. 2003; Nagao et al. 2006a; see 



also Juarez et al. 2009). Table 3 lists the inferred emission-line 
flux ratios. The error on the emission-line flux ratios is estimated 
by taking account both statistical and systematic errors. The sys- 
tematic error (about five percent of the line flux) is associated 
with the uncertainties in the line-profile fitting. 
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Fig. 3. Each panel shows the composite spectra in a fixed interval of M B h by varying the Eddington ratio. The height scale is the 
same for all spectra, but different spectra are shifted by two. 



To investigate the M bh -Z B lr relation and L/L E dd-Z B LR re- 
lation, we show the dependence of the emission-line flux ra- 
tios on M B h and L/L^i in Figure 5. The left-hand panels in 
Figure 5 show that there is a tight correlation between the 
metallicity-sensitive emission-line flux ratios and M B h- On the 
other hand, the right-hand panels show that there is a large scat- 



ter and no clear dependence of the emission-line flux ratios on 
the Eddington ratio. These results are consistent with the earlier 
findings of Warner et al. (2004), who analyzed rest-frame ultra- 
violet spectra of type 1 AGNs and found the correlation between 
M B h and Z B lr- In contrast, these results appear to be inconsis- 
tent with the results of Shemmer et al. (2004), who found a sig- 
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9.6Slog(jl/ BH /M e )<9.B 



9.8glog(M BH /tf o )<10.0 




rest 



(A) 



rest 



(A) 



Fig. 3. Continued. 



nificant relation between L/L E dd and the N v/C iv line flux ratio 
while the dispersion in the relation between M B h and N v/C iv 
is much larger (see also Dietrich et al. 2009). It should be noted 
that the sample of Shemmer et al. (2004) includes AGNs both 
at high-z and in the local universe, and the large dispersion 
in the relation between M B h and Nv/Civ in their sample is 
mainly caused by local narrow-line Seyfert 1 galaxies (NLSls). 
However, it is known that NLSls do not follow usual trends 
in terms of BLR emission-line flux ratios (Shemmer & Netzer 
2002), that is also shown by Warner et al. (2004). Actually, 
the dispersion in the M B H-Nv/Crv relation in the sample of 
Shemmer et al. (2004) is significantly reduced if the NLSls are 
removed. Note that most of the high-MBH objects in the sample 
of Shemmer et al. (2004) are very luminous quasars and with 
high-L/L E dd, because their high-z sample was observed spectro- 
scopically in the near-infrared and therefore high-luminosity ob- 
jects were preferentially selected. Therefore, due to the resulting 
strong degeneracy between Mbh and L/L E dd, the N v/Civ flux 
ratio in the sample of Shemmer et al. (2004) shows significant 
dependence on both M BH and L/L E dd- 



Table 4. Results of the Spearman rank-order test. 



line ratio 



P(r s ) 1 



vs. black hole mass 



N v/C iv 
Nv/Hen 
(Siiv+Oiv])/Civ 
Aim/Civ 



+0.94 
+0.94 
+0.99 
+0.87 



5.03 x 10~ 15 
9.08 x 10~ 15 
2.20 x 10- 16 
3.44 x 10-'° 



vs. Eddington ratio 



N v/C iv 
Nv/Hen 
(Siiv+Oiv])/Civ 
Aim/Civ 



+0.86 
+0.97 
-0.28 
+0.38 



5.09 x 10" 10 
2.20 x 10- 16 
1.24 x 10-' 
3.60 x 10- 2 



" Spearman rank-order correlation coefficient. 

* Probability of the data being consistent with the null hypothesis 

that the flux ratio is not correlated with the black hole mass or the 

Eddington ratio. 



In order to further investigate the dependence of Zblr on 
M B h and L/L E dd, it is crucial to examine the M Bll -Z BLR rela- 
tion and the L/L E dd-Z BE R relation independently . For instance, 
if we plot the Nv/Civ flux ratio as a function of L/LEdd re- 
gardless of M B h, possible correlations between Nv/Civ and 
Mbh would prevent us from investigating the relation between 
N v/C iv and L/L E dd- In order to minimize the effect of L/LEdd 
or Mbh on any relation between the emission-line flux ratios 
and M B h or L/L E dd respectively, we normalize the emission-line 
flux ratios by the values measured from the composite spectra 
at -0.6 < log(L/L Ed d) < -0.4 or 9.2 < log(M BH /M ) < 9.4, 
respectively. In Figure 6, we show the normalized flux ratios as 
a function of Mbh or L/L E dd- In the left-hand panels of Figure 6, 
we find that there are clear relations between M B h and emission- 
line flux ratios as seen in Figure 5. On the other hand, surpris- 
ingly, we also find that there is a correlation between L/L E dd 
and the emission-line flux ratios involving N v in the right pan- 
els, which was not obvious in Figure 5. To assess these relations 
quantitatively, we apply the Spearman rank-order test to statisti- 
cally examine the possible correlations between Zblr and Mbh, 
or L/LEdd- We calculate the Spearman rank-order correlation co- 
efficient r s and their statistical significance p(r s ), which are given 



in Table 4. The results of the Spearman rank-order test are sum- 
marized as follows: 

- there are significant correlations between the black hole mass 
and all four emission-line ratios, and 

- the line flux ratios involving N v (i.e., N v/C iv and N v/He n) 
are positively correlated with the Eddington ratio, while the 
other two flux ratios, i.e., (Si iv+O rv])/C iv and Al m/C iv, do 
not show any significant correlation with L/L E dd- 



5. Discussion 

5.1. The black hole mass-metallicity relation 

As presented in section 4 (see also Figures 5 and 6), the tight 
relation between Mbh and all metallicity sensitive line ratios in- 
dicates that the L A gn-Zblr correlation is driven mostly by the 
Mbh-Zblr relation (and not by any L/L E dd-ZBLR relation). The 
Mbh-Zblr relation can be understood in the following terms. If 
the M B H-Mbui relation observed in the local universe (Magorrian 
et al. 1998) holds at z ~ 2 - 3 (Janke et al. 2009), the M B h-Z B lr 
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Fig. 4. Composite spectra of each quasar subsample and fitting results, shown in the logarithmic scale (upper panels). The residuals 
are shown in the lower panels. Dark grey lines show the wavelength regions used for the fitting, while the wavelength regions 
excluded from the fitting process are shown in light grey. The spectral regions denoted by vertical dashed lines are used to determine 
the initial guess of the continuum level and slope in the fitting process. The initial guess and the fitted power-law continuum are 
denoted by red dotted and dashed lines, respectively. 
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Fig. 4. Continued. 
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Table 2. The line profile and the velocity shift of HILs and LILs. 


Black Hole Mass 


Eddington Ratio 


Line 


Line Profile 


Velocity Shift " 










a 


P 


[km/s] 


8.4 < log(M BH /M G ) < 8.6 


-0.4 < log(L/L Edd ) 


< -0.2 


HIL 


123.02 ± 0.05 


110.88 ±0.04 


-605.6 ± 0.4 








LIL 


138.66 ±0.11 


178.94 ± 0.36 


-281.4 ± 0.5 




-0.2 < log(L/L Edd ) 


< 0.0 


HIL 


165.74 ± 0.06 


123.38 ± 0.04 


64.5 ± 0.5 








LIL 


141.34 ±0.11 


110.46 ±0.18 


-159.2 + 0.5 


8.6 < log(M BH /M ) < 8.8 


-0.6 < loe(L/L F dd) 

— * V ^ \*— ' / *CUU / 


< -0.4 


HIL 


127.87 ± 0.06 


101.77 ± 0.04 


-266.4 ± 0.6 








LIL 


118.05 ±0.10 


140.65 ± 0.29 


-248.8 ± 0.6 




-0.4 < log(L/L Fdd ) 


< -0.2 


HIL 


132.01 ± 0.06 


116.01 ± 0.04 


-455.8 ± 0.3 








LIL 


122.14 ±0.12 


120.84 ± 0.27 


81.5 ± 1.3 




-0.2 < log(L/L Edd ) 


< 0.0 


HIL 


129.78 ± 0.06 


122.62 ± 0.06 


-287.8 ± 0.6 








LIL 


117.38 ±0.12 


109.67 ± 0.22 


-3.5 ± 1.3 


8.8 < log(M BH /M ) < 9.0 


-0.8 < log(L/L Edd ) 


< -0.6 


HIL 


81.28 ±0.05 


101.93 ±0.04 


-444.2 ± 0.8 








LIL 


81.38 ±0.08 


69.78 ±0.11 


-96.5 ± 2.7 




-0.6 < log(L/L Fdd ) 

w — ""cV / cut!/ 


< -0.4 


HIL 


103.79 ± 0.05 


101.96 ±0.05 


-668.9 ± 0.3 








LIL 


119.18 ± 0.11 


123.08 ± 0.34 


227.4 ± 1.2 




-0.4 < log(L/L Fdd ) 


< -0.2 


HIL 


121.21 ±0.06 


104.21 ± 0.05 


-381.9 ±0.7 








LIL 


116.67 ±0.12 


106.54 ± 0.21 


97.4 ± 1.4 




-0.2 < log(L/L Fdd ) 

" — Ov ' CUU/ 


< 0.0 


HIL 


142.79 ± 0.07 


113.59 ±0.05 


-276.5 ± 0.7 








LIL 


124.61 ±0.12 


106.89 ± 0.25 


147.0 ± 1.5 


9.0 < log(M BH /M ) < 9.2 


l.KJ 1^ lUg^L/ l^ Edd J 


< -0.8 


HIL 


77 w + n ()S 




-520.0 ± 0.9 








LIL 


90.39 ±0.10 


49.41 ± 0.07 


988.6 ±0.9 




-0.8 < log(L/L Edd ) 


< -0.6 


HIL 


95.05 ± 0.06 


86.71 ± 0.05 


-470.5 ± 1.0 








LIL 


121.39 ±0.14 


96.99 ±0.19 


548.1 ± 1.3 




-0.6 < log(L/L Fdd ) 


< -0.4 


HIL 


98.13 ±0.06 


92.36 ± 0.05 


-316.5 ± 1.0 








LIL 


102.07 ±0.11 


87.34 ±0.16 


391.4 ± 1.4 




-0.4 < log(L/L Edd ) 


< -0.2 


HIL 


108.65 ± 0.07 


95.04 ± 0.05 


-169.0 ± 1.0 








LIL 


105.07 ±0.12 


93.23 ±0.18 


472.7 ± 1.6 




-0.2 < logfL/Lpdd) 


< 0.0 


HIL 


129.39 ± 0.07 


100.76 ± 0.06 


-315.6 ±0.8 








LIL 


111.97 ±0.12 


120.67 ±0.31 


148.6 ± 1.5 


9.2 < log(M BH /M ) < 9.4 


i.u ^ iog^/.L Edd ,j 


< -0.8 


HIL 


SI AT, + Oh 


OJ.7U U.UU 


-586.6 ± 1.2 








LIL 


98.97 ±0.14 


77.96 ±0.16 


1313.4 ± 1.4 




-0.8 < log(L/L Fdd ) 


< -0.6 


HIL 


88.26 ± 0.07 


75.74 ± 0.05 


-483.7 ± 1.2 








LIL 


88.12 ±0.12 


91.19 ±0.21 


605.9 ± 1.5 




-0.6 < log(L/L Fdd ) 


< -0.4 


HIL 


91.99 ±0.07 


79.27 ± 0.05 


-377.1 ± 1.1 








LIL 


94.80 ±0.12 


86.50 ±0.21 


504.4 ± 1.8 




-0.4 < log(L/L Fdd ) 


< -0.2 


HIL 


101.37 ±0.08 


86.08 ± 0.06 


-326.4 ±1.1 








LIL 


93.55 ± 0.13 


94.11 ± 0.24 


569.3 ± 1.8 




-0.2 < log(L/L Fdd ) 


< 0.0 


HIL 


134.51 ± 0.09 


100.59 ± 0.06 


-36.7 ± 0.7 








LIL 


109.86 ±0.15 


144.88 ± 0.41 


227.6 ± 1.2 


9.4 < log(M BH /M ) < 9.6 


1 - z — Iu gV L / -MidtU 


< -1.0 


HIL 


fS1 10 + 05 


71 61 + 05 


-1066.5 ± 1.7 








LIL 


73.61 ±0.10 


37.91 ±0.04 


218.0 ±2.3 




-1.0 < log(L/L Edd ) 


< -0.8 


HIL 


73.84 ± 0.07 


69.71 ±0.05 


-367.1 ± 1.6 








LIL 


85.60 ±0.13 


63.20 ±0.16 


1040.9 ± 2.0 




-0.8 < log(L/L Edd ) 


< -0.6 


HIL 


72.02 ± 0.07 


66.98 ± 0.06 


-503.1 ± 1.6 








LIL 


84.97 ±0.13 


63.83 ±0.18 


1031.0 ± 2.3 




-0.6 < log(L/L Fdd ) 


< -0.4 


HIL 


95.55 ± 0.09 


71.92 ±0.06 


-195.7 ± 1.4 








LIL 


95.18 ±0.15 


80.44 ± 0.20 


1108.0 ± 1.8 




-0.4 < log(L/L Fdd ) 


< -0.2 


HIL 


101.12 ±0.08 


76.95 ± 0.06 


-354.3 ± 1.3 








LIL 


102.99 ±0.15 


92.29 ± 0.22 


889.8 ± 1.7 


9.6 < log(M BH /M ) < 9.8 


i.z, ^ lug^L/ i^ Edd ; 


< -1.0 


HIL 


f,^ 87 + n Q7 


64 40 + 06 


-975.1 ± 1.9 








LIL 


100.74 ±0.19 


51.78 ±0.17 


1120.5 ±2.2 




-1.0 < log(L/L Fdd ) 

— Ov ' CUU/ 


< -0.8 


HIL 


67.35 ± 0.08 


52.80 ± 0.05 


-467.1 ± 1.8 








LIL 


83.13 ±0.16 


53.26 ±0.16 


1415.7 ± 2.4 




-0.8 < log(L/L Fdd ) 


< -0.6 


HIL 


70.16 ± 0.08 


58.12 ± 0.06 


-491.7 ± 1.8 








LIL 


83.86 ±0.16 


55.64 ± 0.16 


1384.7 ± 2.4 




-0.6 < log(L/L Edd ) 


< -0.4 


HIL 


81.89 ±0.08 


64.19 ±0.06 


-538.3 ± 1.7 








LIL 


93.55 ±0.16 


68.60 ±0.19 


1184.9 ± 1.8 




-0.4 < log(L/L Edd ) 


< -0.2 


HIL 


92.17 ±0.08 


70.41 ± 0.06 


-528.2 ± 1.3 








LIL 


94.59 ±0.14 


86.38 ± 0.23 


773.8 ± 1.8 


9.8 < log(M B H/M ) < 10.0 


-1.2 < log(L/L Edd ) 


< -1.0 


HIL 


60.67 ± 0.09 


44.19 ±0.06 


nor\ ft i t ft 
680.0 ± 1.9 








LIL 


75.89 ±0.16 


48.30 ±0.14 


1723.7 ± 2.2 




-1.0<log(L/L Edd ) 


< -0.8 


HIL 


65.97 ± 0.09 


48.31 ±0.06 


-396.6 ±2.1 








LIL 


73.31 ±0.17 


43.66 ±0.12 


1995.4 ± 2.2 




-0.8 < log(L/L Edd ) 


< -0.6 


HIL 


73.09 ± 0.09 


47.22 ± 0.06 


-208.3 ± 1.9 








LIL 


79.10 ±0.18 


46.35 ±0.15 


2018.2 ± 2.6 




-0.6 < log(L/L Edd ) 


< -0.4 


HIL 


82.84 ± 0.09 


51.92 ±0.06 


-585.6 ± 1.6 








LIL 


97.84 ±0.19 


63.54 ± 0.20 


1480.2 ± 1.9 



" The shift of the emission line center (A ) from the systemic velocity. 
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relation obtained in this study implies a M\, u \-Zblr relation. In 
addition, if the BLR metallicity is associated with the metal- 
licity of the host galaxy, the Mbh-Zblr relation would simply 
trace the well-known correlation between the mass and metallic- 
ity in galaxies, i.e., the Mb u i-Zb u i relation (Lequeux et al. 1979; 
Tremonti et al. 2004; Maiolino et al. 2008). This implies that 
the origin of the Lagn-Zblr relation is the M\, u \-Zb u \ relation (or, 
more simply, the "mass-metallicity relation"). 

The one potential problem of this interpretation is the lack of 
redshift evolution of the Lagn-Zblr, out to z ~ 6 (Nagao et al. 
2006a; Juarez et al. 2009; see also Dietrich et al. 2003), which 
is at odds with the redshift evolution of the mass-metallicity 
relation observed in galaxies (Maiolino et al. 2008, Mannucci 
et al. 2009). There are two possible solution to this problem. 
More recently, Mannucci et al. (2010) have found that the mass- 
metallicity relation in galaxies also depends on the SFR, more 
active galaxies being characterized by lower metallicity. They 
also found that, such mass-metallicity-SFR relation does not 
evolve with redshift out to z ~ 2.5; previous claims on the evolu- 
tion of the mass-metallicity relation were simply a consequence 
of the fact that high-z surveys have probed galaxies on average 
with higher SFR at high redshift. If quasar host galaxies at differ- 
ent redshifts are characterized by similar SFR (e.g., Maiolino et 
al. 2007a; Netzer et al. 2007; Lutz et al. 2008), this would imply 
that they are likely characterized by the mass-metallicity rela- 
tion. However, this explanation is not totally satisfactory, since 
actually there is evidence that even the mass-metallicity-SFR re- 
lation evolves at z > 3 (Mannucci et al. 2010), while the Lagn- 
Zblr in quasars does not evolve out to z ~ 6. 

Another effect, which goes in the direction of making the 
non-evolution of the Lagn-Zblr relation even more problematic, 
is that at high redshift the MBH-A^bui relation seems to evolve to- 
wards higher Mbh/M\, u \ relative to the local relation (Walter et 
al. 2004; Woo et al. 2006, 2008; Maiolino et al. 2007b; Lamastra 
et al. 2010; Merloni et al. 2010). As a consequence, for a given 
black hole mass (hence similar luminosities, for a given accre- 
tion rate) high redshift AGNs should be associated with less 
massive, hence lower metallicity, host galaxies relative to AGNs 
at lower redshift. 

A possible redshift evolution of the accretion rate L/Lndd (at 
a given black hole mass, Netzer & Trakhtenbrot 2007), would 
make the problem of the non evolution of the Lagn-Zblr rela- 
tion even more serious. Indeed, at a given luminosity a high-z 
AGN would have a black hole mass lower than a lower redshift 
AGN, hence should be associated with a less massive, and less 
enriched, host galaxy. 



A possible explanation is that luminous quasars host galaxies 
do not follow the same secular evolutionary process as most of 
the galaxies used to investigate the metallicity evolution, but are 
instead characterized by vigorous star formation likely associ- 
ated with merging events (e.g., Shao et al. 2010), yielding a very 
rapid co-evolution of the black hole and its host galaxy, imply- 
ing a very rapid metals enrichment. Within such strong an rapid 
co-evolutionary scenario it is likely that the metallicity investiga- 
tion of quasars host galaxies through the cosmic epochs is sub- 
ject to strong selection effects. More specifically, as discussed 
in detail in Juarez et al. (2009), because of the co-evolutionary 
link between black holes and their host galaxies, quasars are de- 
tectable in flux-limited surveys only once their host galaxies are 
already evolved and chemically enriched, and this connection is 
the same at any epoch, regardless of redshift (see also Kawakatu 
etal. 2003). 

5.2. The dependence between the Eddington ratio and 
emission-line ratios involving N v 

In section 4, we have also found that the emission-line flux ra- 
tios involving Nv are related to L/Lsdd whereas the other two 
emission-line flux ratios, i.e., (Siiv+Oiv])/Civ and Aim/Civ, 
do not show any dependence on L/Lndd- If the relation between 
the Eddington ratio and emission-line flux ratios involving N v is 
ascribed to a L/Ledd-ZfiLR relation, then the other two emission- 
line flux ratios should also show dependence on L/Lndd, since all 
four emission-line flux ratios are sensitive to the gas metallicity 
of the BLR. These results thus indicate that emission-line flux ra- 
tios involving N v do not simply depend on metallicity but also 
on other independent parameters. One possible idea is that the 
difference of the behavior between flux ratios with and without 
N v may be associated with the different enrichment timescale of 
nitrogen relative to other elements. 

Recent studies on local AGNs have shown that nuclear star 
formation and AGN activity are not coeval (Davies et al. 2007; 
see also Wild et al. 2010). More specifically, black hole accre- 
tion appears to occur in post-starburst nuclei, with a delay of 
~ 10 s years from the starburst. Probably during the active star- 
burst phase powerful supernova explosions expel the circumnu- 
clear gas preventing it to reach the black hole, while the more 
gentle winds of AGB stars occurring on timescales of 10 8 years 
are capable of stirring the interstellar matter (ISM), making it 
loose angular momentum and then feed the AGN. The produc- 
tion of nitrogen is delayed with respect to other elements, and 
indeed produced mostly by AGB stars whereas other elements 
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Table 3. The measured line flux ratios. 



Rlack Hole Mass 




Eddington Ratio 




N" v/P TV 


Line Flux Ratios " 
Nv/Hen (Siiv+Oiv])/Civ 


Al TTT/P TV 

r\i 111/ \^ IV 


8.4 < log(M BH /M Q ) < 


3.6 


-0.4 < log(L/L Edd ) 


< 


-0.2 


0.320 ± 0.024 


1.974 ±0.144 


0.197 ±0.015 


0.049 ± 0.004 






-0.2 < log(L/L Edd ) 


< 


0.0 


0.460 ± 0.034 


3.556 ±0.259 


0.193 ±0.013 


0.045 ± 0.003 


8.6 < log(M BH /M Q ) < 


3.8 


-0.6 < loe(L/L FH ,i) 


< 


-0.4 


0.489 ± 0.037 


2.900 ± 0.212 


0.237 ±0.018 


0.055 ± 0.004 






-0.4 < log(L/L Edd ) 


< 


-0.2 


0.480 ± 0.035 


3.170 ±0.231 


0.214 ±0.014 


0.049 ± 0.004 
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-0.8 < log(L/L Edd ) 
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-0.6 < log(L/L Edd ) 
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0.118 ±0.010 



" Line flux ratios of Nv^1240/Civ^1549, Nv^l240/He 1^1640, (Siiv^l398+Oiv]^1402)/Civil549, and A1iil11857/Civ,11549. These errors 

are estimated by taking into account both random error and systematic error (five percent on line flux) for each emission line. 
* In this interval Al ni/11857 is not detected. 



(i.e., Si, O, and Al) are produced by massive stars on much 
shorter timescale. Carbon originates mainly from low-mass stars 
(Chiappini et al. 2003); their life time is ~ 10 9 ~ 10 years. This 
means that carbon would not be associated with black hole ac- 
cretion. Moreover, since Civ/11549 is a strong coolant of the 
BLR, its intensity is little sensitive to the carbon abundance. As 
a consequence, if also at high-z black hole accretion is associ- 
ated with the delayed onset of AGB stars, one naturally expects 
a correlation between Eddington ratio and nitrogen abundance, 
as observed by us. Note that, in this picture, the correlation be- 
tween Eddington ratio and emission-line flux ratios involving 
N v suggest that the timescale of AGN feeding are shorter than 
the timescale of the nitrogen enrichment by AGB stars (< 10 8 
years). 

The relation between L/LEdd and nitrogen abundance is as- 
sociated to black hole accretion and star formation, i.e. on phe- 
nomena occurring on relatively short time scales (~ 10 8 years). 
Conversely, the relation between black hole mass and metallicity 
is more fundamental and representative of the global evolution 



of these systems, since it connects physical quantities integrated 
over the whole formation history of black holes and galaxies. 

Within this context it is interesting to note that NLSls, 
which are characterized by very high Eddington ratios, have high 
N v/C iv ratio (Shemmer et al. 2004), in line with the trend found 
above for luminous high-z quasars and suggesting that NLSls 
have undergone vigorous star formation (Nagao et al. 2002). 

It should be noted that previous studies have found that 
quasar hosts are characterized by vigorous ongoing star forma- 
tion (e.g., Maiolino et al. 2007a; Netzer et al. 2007; Lutz et al. 
2008). Our finding of a correlation between black hole accre- 
tion rate and nitrogen abundance suggests that at earlier times, 
by about 10 8 years, star formation was even higher. 



6. Conclusions 

We have investigated the relationship between the metallicity 
and SMBH mass or Eddington ratio by producing composite 
SDSS quasar spectra with 33 subsamples divided in intervals of 
black hole mass and Eddington ratio at 2.3 < z < 3.0. In each of 
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Fig. 5. Relation between the metallicity sensitive emission-line flux ratios and the black hole mass (left hand panels), and the 
Eddington ratio (right hand panels). 



these spectra we measure emission-line flux ratios that are sen- 
sitive to the metallicity of the BLR gas. We have investigated the 
Mbh-Zblr dependence or the L/Lem-Zblr dependence by com- 
paring the emission-line flux ratios of the composite spectra by 
fixing the Eddington ratio or black hole mass, respectively. We 
have found the following results: 



1 . There is a tight correlation between the BLR metallicity and 
black hole mass, while the BLR metallicity is almost in- 
dependent of the Eddington ratio, implying that the Lagn- 
Zblr relation in quasars results from the relationship be- 
tween Zblr and Mbh (and not from a relationship between 
Z BL r and L/L Edd ). 
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Fig. 6. Metallicity sensitive emission line flux ratios normalized 
by the value measured in the composite spectra with -0.6 < 
log(L/LEdd) < -0.4 for the individual black hole mass bins, as a 
function of black hole mass (left hand panels). Emission line flux 
ratios normalized by the value at 9.2 < log(MBH/A^o) < 9.4 for 
the individual Eddington ratio bins, as a function of Eddington 
ratio (right hand panels). Symbols are the same as those in Figure 
5. 



2. The emission-line flux ratios involving Nv (i.e., Nv/Civ 
and N v/He n) correlate with the Eddington ratio (while the 
other two metallicity sensitive emission-line flux ratios do 
not show any correlation with L/Lndd)- This result sug- 
gests that the accretion rate in quasars is associated with a 
post-starburst phase, when AGB stars enrich the interstellar 
medium with nitrogen. 

By assuming that the Mbh-MmI relation applies at high-redshift, 
the Mbh-Zblr relation is likely a consequence of the M\, u \-Zb u \ 
relation in galaxies. 



The lack of redshift evolution of the Lagn-Zblr relation does 
not necessarily imply a lack of metallicity evolution in AGN host 
galaxies. Indeed, as already suggested in previous works (e.g., 
Juarez et al. 2009), a combination of selection effects and of the 
co-evolution between black hole and galaxies causes quasars to 
be detected only once their host galaxy is already chemically 
evolved at any epoch, regardless of redshift. 
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